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ABSTRACT 

Previous  considerations  by  asymptotic  expansion  procedures  of  the 
relation  between  elasticity  theory  results  and  thin-shell  theory  results  for 
the  case  of  rotationally  symmetric  deformations  of  an  edge -loaded  semi- 
infinite  circular  cylindrical  shell  are  supplemented  by  an  analysis  of  this 
problem  for  a shell  possessing  a limiting-type  orthotropy,  such  that  transverse 
normal  strains  vanish  identically.  It  is  shown  that  assuming  this  kind  of 
orthotropy  has  the  important  benefit  of  allowing  the  derivation  of  exact  expres- 
sions for  the  edge  zone  solution  contribution,  when  such  exact  expressions 
are  not  possible  for  the  problem  of  the  shell  with  more  general  properties 
of  the  material.  — One  result  of  our  analysis  is  an  answer  to  the  following 
question.  Given  a shell  with  arbitrarily  prescribed  edge  displacements,  what 
is  the  asymptotically  exact  form  of  the  corresponding  conditions  for  this  same 
problem,  treated  within  the  framework  of  two-dimensional  thin-shell  theory? 
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ON  TWO-  AND  THREE-DIMENSIONAL  RESULTS 
FOR  ROTATIONALLY  SYMMETRIC  DEFORMATIONS 
OF  CIRCULAR  CYLINDRICAL  SHELLS 

S.  Nair  and  E.  Reissner 


INTRODUCTION 


We  return  once  more  to  the  problem  of  rotationally  symmetric 
deformations  of  a semi-infinite  circular  cylindrical  shell  as  the  simplest 
non-trivial  example  of  the  relation  between  three-dimensional  elasticity- 
theory  analysis  and  two-dimensional  thin  shell-theory  analysis.  The  first 
paper  on  this  subject  [l]  considered  the  problem  of  the  asymptotic  determina- 
tion of  a class  of  "interior"  solutions  for  the  given  three-dimensional  boundary 
value  problem  and  the  derivation  therefrom  of  a system  of  two-dimensional 
shell-theory  equations  including  the  formulation  of  shell-theory  boundary 
condition  statements  from  given  three-dimensional  elasticity-theory  state- 
ments cf  such  conditions,  for  the  case  that  these  conditions  were  stress 
boundary  conditions.  A subsequent  paper  by  Reiss  [Z]  extended  this  work  by 
considering  complete  asymptotic  solutions,  including  interior  solutions  and 
"edge-zone"  solution  contributions.  The  results  obtained  in  this  manner 
confirmed  the  conclusions  in  [l]  in  regard  to  the  problem  of  two-dimensional 
shell  theory,  while  at  the  same  time  supplying  significant  additional  insights 
in  regard  to  the  nature  of  the  relation  between  two-  and  three-dimensional 
theory,  with  these  insights  having  meanwhile  been  extended  and  generalized 
in  important  ways  by  various  other  workers,  in  particular  by  Goldenweiser  [3]. 


One  of  the  difficulties  encountered  in  the  use  of  an  edge -zone  solution 

contribution,  as  done  by  Reiss  [2],  consists  in  the  fact  that  the  relevant  two- 

dimensional  boundary  value  problem  for  a bi-harmonic  function  defined  in  a 

semi-infinite  strip  cannot,  for  some  important  cases  including  the  case  of  pure  fion 

on  □ 

traction  conditions  and  of  pure  displacement  conditions  be  solved  in  closed  form,  n 


and  to  the  extent  that  this  is  the  case  the  asymptotic  results  which  are  obtained 
remain  approximate  rather  than  exact. 
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Given  the  impossibility  of  a closed-form  solution  of  the  relevant 
bi-harmonic  problem,  as  well  as  the  apparent  absence  of  results  for  the 
case  of  pure  displacement  edge  condition  cases,  we  have  recently 
considered  the  complete  problem  by  combining  interior  asymptotic 
expansions,  Rayleigh-Ritz  type  edge-zone  solution  contributions,  and 
upper  and  lower  bound  formulas  through  use  of  the  principles  of  minimum 
potential  and  complementary  energies  [4].  The  principal  results  of  this 
analysis  consisted  in  the  derivation  of  upper  and  lower  bounds  for  the  values 
of  influence  coefficients  involved  in  the  solution  of  the  semi-infinite  circular 
cylindrical  shell  problem  with  prescribed  edge  tractions  or  prescribed  edge 
displacements.  All  of  these  bound  results  were  such  as  to  imply  the  deter- 
mination of  exact  results  for  the  solution  of  the  three-dimensional  problem 
by  means  of  two-dimensional  theory,  insofar  as  the  leading  terms  in  an 
expansion  of  the  solution  of  the  three-dimensional  problem  in  powers  of  wall 
thickness  h to  shell  radius  ratio  a were  concerned.  We  also  determined 
supplementary  terms  for  such  an  expansion,  of  relative  order  h/a  (including 
terms  of  order  (h/a)^  which  are  encountered  for  some  classes  of  edge  condi- 
tions), with  these  supplementary  bound  terms  being  such  that  in  some  cases 
there  was  coincidence  between  upper  and  lower  bound  results  so  that,  in 
effect,  a determination  of  exact  results,  up  to  terms  of  relative  order 
h/a,  was  accomplished. 

In  order  to  understand  the  meaning  of  these  results,  as  well  as  the 
significance  of  the  analysis  which  follows  it  is  convenient  to  interpret  the 
leading-term  results  as  exact  results  for  an  infinitely-thin  shell  (that  is,  for 
a shell  for  which  h/a-»0),  with  the  supplementary  terms  representing  the 
effects  of  finite  thickness.  There  are  altogether  three  distinct  effects  of 
finite  thickness.  The  first  of  these  is  a geometrical  effect,  having  to  do  with 
the  change  of  width  of  shell  elements  with  distance  from  the  middle  surface 
(so  that  this  effect  is  absent  for  the  special  case  of  a flat  plate).  It  is  known 
that  this  effect  is  taken  account  of  properly  (assuming  absence  of  the  other 
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two)  in  a refined  two-dimensional  shell  theory  associated  with  the  names  of 
Flvigge,  Lurie  and  Byrne.  The  second  effect  is  the  effect  of  transverse 
normal  stress  deformability  (vhich  is  absent  for  the  case  of  a limiting-type 
orthotropic  material  unable  to  sustain  transverse  normal  strains).  It  has 
earlier  been  shown  that  this  effect  is  of  the  same  order  of  magnitude  as  the 
geometrical  FLB  effect  [l],  and  our  recent  work  [4]  indicates  that  for  some 

i 

classes  of  edge  condition  this  effect  comes  out  to  be  of  relative  order  (h/a)^ 
(without  these  terms  being  of  numerical  significance,  however,  in  comparison 
with  the  co-existing  h/a-order  effects). 

The  third  of  the  effects  of  finite  thickness  is  the  effect  of  transverse 
shear  deformability.  This  effect  too  comes  out  to  be  of  relative  order  h/a. 

We  are  not  concerned  here  with  the  consequences  of  this  effect,  in  regard  to 
the  order  of  the  differential-equation  system  and  to  the  number  of  the  asso- 
ciated boundary  conditions,  as  discussed  most  simply  in  recent  work  dealing 
with  the  subject  of  plates  [5].  Rather,  we  are  concerned  with  this  effect 
from  the  point  of  view  of  its  relative  numerical  dominance  in  comoarison 
with  the  other  two,  as  revealed  by  our  upper  and  lower  bound  calculations  [4]. 

Having  previously  obtained  upper  and  lower  bound  results  for  the 
three  effects  of  geometry,  transverse  normal  stress  deformability,  and 
transverse  shear  deformability,  with  these  three  effects  being  additive  up  to 
orders  of  magnitude  which  are  of  primary  interest,  we  now  undertake  an 
asymptotic  analysis  of  two  of  the  three  effects,  these  being  the  effects  of 
transverse  shear  deformability  and  of  geometry.  Our  analysis  is  based  on 
recognition  of  the  fact  that  it  is  possible  to  derive  exact  solutions  for  the 
two-dimensional  semi-infinite  strip  problem  governing  the  edge-zone 
solution  contribution,  upon  assuming  a limiting -type  orthotropy  in  such  a 
manner  that  transverse  normal  strains  vanish  identically. 

Having  the  existence  of  these  exact  solutions  for  the  edge-zone 
contributions  involved  in  the  asymptotic  expansion  procedure,  we  are 
now  in  a position  to  verify  and,  in  principle,  to  refine  the  results 
of  our  upper  and  lower  bound  analysis.  Beyond  this,  we  are  able  to  obtain 
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results  for  types  of  boundary  conditions  which  do  not  fall  within  the  scope 
of  the  indicated  bound  solutions.  To  mention  a specific  example,  our  analysis 
permits  us  to  solve  a problem  which  has  long  been  of  interest  to  us  but  for 
which  until  now  no  rational  solution  has  come  to  our  attention.  The  problem 
is  as  follows.  Given  a semi-infinite  circular  cylindrical  shell,  with  arbitrarily 
prescribed  edge  displacements  as  loading  conditions.  To  be  determined  is  the 
asymptotically  exact  form  of  the  corresponding  conditions  of  the  first-order 
interior  solution  contribution,  to  wit  the  appropriate  form  of  the  corresponding 
boundary  conditions  for  this  same  problem,  treated  within  the  frame  work  of 
standard  two-dimensional  thin  shell  theory. 

FORMULATION  OF  THE  PROBLEM 


We  take  as  differential  equation  for  symmetrical  deformations  of 
circular  cylindrical  bodies  a system  consisting  of  the  equilibrium  equations 

"x,x ° = ° • (2-1) 

in  conjunction  with  stress  strain  (displacement)  relations  of  the  form 

<y  <7^  - I/O  O 
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u = - V — — 
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u + V = — , 

, r , X G ’ 


(2.2) 


where  E^  = (EE^)^,  with  positive  E,  E^  and  G,  and  with  the  additional  strain 
energy  positive -definiteness  conditions  P®  < 1 and  2p®  < 1 - p. 

The  system  (2.  1)  and  (2.  2)  is  to  be  solved  in  the  region  0 ^ x < ®, 
a-cSffia+c  subject  to  face  boundary  conditions 


a ± c; 


a =0 
r 


T = 0 , 


(2.3) 


subject  to  edge  boundary  conditions  "at  infinity"  which  for  all  cases  are  taken 
in  the  form 


X -♦  “S 


u = 0 . 


V = 0 


(2.4) 
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and  subject  to  edge  boundary  conditions  at  the  loaded  edge  of  the  shell,  of 
the  form 


(2.5) 


with  the  right  hand  sides  in  these  four  relations  being  presc ribed  functions 

of  r,  subject  only  to  the  restriction  that /^^^(r/aict  dr  = 0. 

a-c  X 

Within  the  frame  work  of  the  above  class  of  three-dimensional  problems 
(which  because  of  the  assumed  absence  of  any  0-dependence  of  the  solutions 
formally  reduce  to  a class  of  two-dimensional  problems)  we  are  particularly 
concerned  in  asymptotic  reductions  to  two-dimensionality  (with  this  reduction 
here  formally  to  one-dimensicnality)  for  the  determination  of  the  weighted 
stress  averages 

a+c  a+c 

M = j (r/a)0’  (r  - a)dr  , Q = J (r/a)Tdr  (2.6) 

a-c  a-c 


and  for  the  determination  of  displacement  measures  such  as, 

V = v(x,  0)  , 0 = u (X,  0)  , (2.7) 

with  these  reductions  being  of  technical  significance  for  "sufficiently  small" 
values  of  the  wall  thickness-diameter  ratio  c/a  of  the  shell. 

In  association  with  the  derivation  of  a system  of  two  (here  one)  -dimen- 
sional differential  equations  for  the  quantities  Q.  V,  0 it  is  necessary  to 

also  derive  a system  of  suitable  edge  conditions,  involving  the  functions 

u » V,  a , T which  appear  in  equations  (2.5).  One  expects,  and  this  has  been 

, r X 

shown  to  be  true  in  [l]  and  [2],  that  insofar  as  the  traction  condition  portions  in 
(2.5)  are  concerned  the  equivalent  lower-dimensional  edge  conditions  are  of  the  form 

X = 0:  M = M , Q = Q (2.8a) 

XX 

with  and  Q defined  in  terms  of  5^  and  T in  accordance  with  equations  (2.6). 

One  also  expects,  and  this  is  generally  considered  to  be  correct,  that 
when  the  displacement  condition  portions  in  equations  (2.  5)  are  such  that  u 
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aa  well  as  v are  independent  of  r,  then  the  equivalent  lower-dimensional 
edge  conditions  are  of  the  form 


X = 0: 


V = V , 


P = u 


(2.8b) 


We  will,  in  what  follows,  re -confirm  the  asymptotic  validity  of  the  above 

i 

expectation,  within  the  frame  work  of  the  restrictions  associated  with  the 
nature  of  the  analysis  which  is  here  carried  out.  Over  and  above  this  result, 
however,  we  will  establish  the  form  of  the  edge  conditions  for  "effective" 
measures  V and  3,  not  necessarily  identical  with  these  measures  as  defined 
in  equations  (2.7),  which  are  valid  in  the  event  that  u and  v are  other  than  indepen- 
dent  of  r.  An  example  of  this  nature  which  will  be  considered  explicitly  is 
the  case  where  u is  assumed  to  be  proportional  to  (r  - a)  in  place  of  the  usually 
assumed  linear  distribution. 

DIFFERENTIAL  EQUATIONS  AND  BOUNDARY  CONDITIONS 
FOR  TRANSVERSELY  INEXTENSIONAL  DEFORMATIONS 

The  case  of  transversely  inextensional  deformations  is  given  upon 

setting 

^r  = “ (3.1) 

in  the  stress  strain  relations  (2.2).  These  may  therewith  be  written  in  the  form 


Er  V = a ^ . VC 
V X 


G(u  + V ) = T , 

$ i.  f A 


(3.2) 


where  it  will  be  assumed  from  here  on  that  E,  V and  G are  independent  of  x and  r. 

Our  first  conclusion  is  now  that  the  transverse  displacement  component 
does  not  vary  across  the  thickness. 


V = V(x)  (3,3) 

and  that,  associated  therewith,  the  transverse  normal  stress  component  c 
assumes  the  character  of  a reactive  quantity. 

In  order  to  solve  the  remaining  boundary  value  problem,  we  begin  by 
satisfying  the  first  of  the  equilibrium  equations  (2.  1)  by  means  of  a stress 
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function 't,  in  the  form 


ra  = , rT  = -a^  , 

X , r ,x 


(3.4) 


where  the  factor  a on  the  right  hand  side  has  been  introduced  to  make  the  writing 
of  some  of  the  developments  which  follow  somewhat  more  convenient. 


We  next  use  the  second  of  the  stress  strain  relations  (3,2),  in  conjunction 
with  equation  (3.3)  in  order  to  write  for  the  circumferential  normal  stress, 

(3.5) 


ra^j  r + EV  . 

0 . r 


With  the  above  expressions  for  (Jq  and  T,  we  obtain  from  the  second 

equilibrium  equation  in  (2,  1)  as  an  expression  for  the  transverse  normal  stress 

CT  -which  satisfies  the  condition  of  vanishing  o for  r = a - c, 
r r 

,r  r 

= a (^  -I-  l/r  )dr  + EVj  r ^dr  . (3.6) 

r J , XX  , r 0 


ra 


a-c 


a-c 


The  condition  that  a^  must  also  vanish  for  r = a4c  gives  as  one  of  two  equations 
connecting  the  two  functions  'I'(x,  r)  and  V(x) 
a 

J ,xx  r , r 

a • c 


a a-c 


= 0 . 


(3.7) 


A second  equation  connecting  and  V follows  upon  expressing  u in  terms 


of  'I'  and  V in  accordance  with  the  last  relation  in  (3.2),  in  conjunction  with  (3.3) 
and  (3.4),  as 


u = -V  — — , 

, r ,x  r G 


(3.8) 


and  by  combining  this  result  with  the  first  relation  in  (3.2),  written  in  the  form 
Eu  = a 


- , with  a and  a«  taken  from  equations  (3.4)  and  (3.  5).  The 

ensuing  differential  equation  may  be  written  in  the  form 


t^v  .iLv  . 0 . 

E \ . rr  r / G , XX  a , XX  ra 


(3.9) 


Having  equations  (3.7)  and  (3.9)  it  remains  to  state  boundary  conditions  in 
terms  of  'I' and  V,  in  accordance  with  the  remaining  face  boundary  conditions 
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t- 


f 

I 


t 


in  (2.  3),  which  concern  T,  and  in  accordance  with  the  edge  conditions  as  stated 
in  Eqs.  (2. 4)  and  (2,  5). 

We  begin  by  satisfying  the  edge  condition  at  infinity  by  stipulating 

'!'(«,  r)  = 0 , V(«)  = 0 . (3.  10) 

With  this,  and  with  observation  of  the  relation  j*  ^^(r/a)C7  dr  = 0,  the  face 

& * c ^ 

boundary  conditions  T(x,  a ± c)  =0  take  on  the  form 

a ± c)  = 0 . (3.11) 

Finally,  the  conditions  (2,  5)  for  the  loaded  edge  of  the  shell  become 


(y  + ~ _ _■ 

or 

>1^ 

, r 

a 

X 

l,x^G  r ” '^,r 

r 

a 

X 

II 

o 

•• 

/ V = V 

or 

a+c 

dr  = -Q 

‘‘a-c  ' 

r 

(3. 12) 


We  note  that  in  writing  the  second  set  of  conditions  in  (3.  12),  we  have  taken 
account  of  the  fact  that  the  assumed  properties  of  the  material  require  that  the 
prescribed  transverse  displacement  v be  independent  of  r and  so  may  be 
written,  in  consistent  fashion,  as  V,  and  that  at  the  same  time  this  specialization 
of  properties  implies  a sensitivity  of  the  medium  to  the  resultant  Q of  the  edge 
stresses  T only,  rather  than  to  the  details  of  a prescribed  T-distribution. 

STIFFNESS  AND  FLEXIBILITY  COEFFICIENTS 


In  accordance  with  our  earlier  work  in  [4],  we  define  general  stiffness 

coefficients  K,  with  reference  to  the  problem  of  prescribed  edge  displacements , 

for  the  case  that  u = -3  and  V = -V  , oy  means  of  relations  of  the  form 
, r o o 

M,(0)  = , Q(0)  = . (4.  1) 

wh.re  Kyj  . Kjy. 

At  the  same  time  we  define  a specialized  system  of  such  coefficients  for 
mixed  edge  condition  problems,  as  follows 
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F 


M^(0)  = (when  T = 0)  , 

Q{0)  = (whena^  = 0)  . 


(4.2) 


► 


Insofar  . :•  the  problem  of  prescribed  edge  tractions  is  concerned,  we 
have  earlier  [4]  defined  flexibility  coefficients  C with  reference  to  traction 
distributions  (r/a)a  = (3M  /2c^)(r  - a)  and  (r/a)T  = (3Q  /4c)[l  - (r  - a)®/c^], 

X O o 

and  weighted  edge  displacement  averages 


= -4^1  t' 

a-c 


^ = -i^J 


,a+c 


[1 


(r  - a)^ 


]vdr 


(4.3) 


a-c 


in  the  form 


P*  = +C..^Q 

o MM  o MQ  o 


V = C^,^M  + C_  _Q 
o QM  o QQ  o 


(4.4) 


with  these  relations  being  "nearly"  the  inverses  of  equations  (4.  1),  upon  the 

identifications  3 =3,  V =V,M(0)=M  and  Q(0)  = Q . 

O O O O X o o 

The  principal  aim  in  [4]  was  the  deduction  of  upper  and  lower  bounds  for  the 

coefficients  K and  C , with  these  bounds  depending  on  the  geometrical  parameter 

c/a,  as  well  as  on  the  materials  property  parameters  E,  E/C,  E/E^,  v and  i> 

Among  the  characteristics  of  these  bound  relations  the  following  are  noteworthy. 

1.  In  the  limit  c/a-»0,  our  upper  and  lower  bound  values  coincide  with 
each  other  and  with  the  corresponding  values  of  what  may  be  called  first-approx- 
imation classical  thin  shell  theory  results. 

2.  For  sufficiently  small  finite  values  of  c/a  our  bound  values 

1 

involve  additive  terms  of  relative  order  (c/a)^  and  of  relative  order  c/a. 

1 

The  terms  of  order  (c/a)^  represent  transverse  normal  strain  effects  alone,  while 
the  terms  of  order  c/a  represent  transverse  normal  strain  effects,  as  well  as 
transverse  shearing  strain  effects,  and  also  geometrical  effects  of  the  Flugge- 
Lurie-Byrne  type.  The  numerically  dominant  of  these  three  effects  is  the  trans- 
verse shearing  strain  effect,  even  though  this  effect  does  not  contain  any 
(c/a)^  - contributions . 


Insofar  as  the  determination  of  stiffness  and  flexibility  coefficients 
is  concerned  the  present  analysis  enables  us  to  obtain  exact  values  of  all 
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[ 


contributions  c5  relative  orders  (c/a)^  and  c/a  (as  well  as  exact  values  of  contri- 
butions of  higher  order  in  c/a  should  we  so  desire)  subject  only  to  the  restriction 
= «.  While  it  is  possible  to  do  this  for  the  entire  set  of  ten  coefficients 
defined  in  equations  (4.  1)  to  (4.4),  the  necessary  analysis  will  be  carried 

out  in  what  follows  only  to  the  extent  of  obtaining  expressions  for  Kqs,  K 

P P Vp 

and  K^. 

THE  NON-DIMENSIONALIZED  BOUNDARY  VALUE  PROBLEM 

We  introduce  a non-dimensional  axial  coordinate  § and  a nondimens ional 
radial  coordinate  7J,  measured  from  the  middle  surface  of  the  shell,  by  writing 

X = b4  , r = a + cT)  , (5.  1) 

with  b being  in  the  nature  of  a characteristic  length,  to  be  chosen  presently. 

We  further  set 

(l-y^)’5'  = T])  , V = V^F(4)  , (5.2) 

where  'I'  and  V remain  to  be  chosen  and  we  define  dimensionless  parameters 
o o ^ 

p and  O'  in  the  form 

P = c/a  , (1  - l/®)ar*  = E/G  . (5.3) 

Indicating  now  differentiations  with  respect  to  ^ and  7}  by  primes  and 
dots,  respectively,  the  two  differential  equations  (3.7)  and  (3.9)  for  'I' and  V 
take  on  the  following  form 


» (g-  — 2 — 

1 + prj 


g*  + + V_[~(l  + PT7)F"  + 


I'P 


o'bs 


1 + prj 


F]  = 0 


(5.4) 

(5.5) 


At  the  same  time  the  face  boundary  conditions  (3.  11)  become 
g(4,  ± 1)  = 0 , (5. 6) 

and  the  edge  boundary  conditions  at  infinity  become 

g(®.  TJ)  = 0 , F(«)  = 0 . (5.7) 
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Insofar  as  the  boundary  conditions  (3.  12)  at  the  loaded  edge  are 
concerned,  we  will  limit  ourselves  in  what  follows  to  the  consideration  of 
just  two  of  a total  of  eight  possible  cases,  both  of  them  concerning  problems 
with  prescribed  axial  displacements  as  expressed  in  terms  of  the  slope  func- 
tion u , with  the  other  condition  being  the  displacement  condition  of  vanishing 
V,  or  the  traction  condition  of  vanishing  Q.  Expressed  in  terms  of  g and  F this 
set  of  conditions  takes  on  the  form 


-u 


. r ' 


pi 

F(0)  = 0 or  J g'(0,  T7)dTj  = 0 


(5.8) 


(5.9) 


Expressions  for  the  stresses  a^,  r,  u q,  and  for  the  stress  resultant  Q and 


the  stress  couple  follow  from  equations  (3.4),  (3.5)  and  (2.6)  as 


<b 


E _o  g* 

“^x  ■ 1 - c 1 + PTJ 


o 


1 - b 1 + PT)  ’ 


V 'I' 

P_£  F E _o  t/g- 

'^8  ■ ^ a 1 + pr?  1 - c 1 + pT)  ' 

E'l'  c _ 'I'  c .1 

o , ^ E o 

= -mjiJ  8'3’j  • Q = — 

-1  -1 


g dl? 


(5.  10) 


(5.  11) 


(5.  12) 


We  list  additionally  as  expression  for  the  circumferential  stress  resultant 
ra+c. 


TT-tr  ^ I+Ott  ^^o^  !‘^  t'gdT? 

^0  ^ FV^p^n  + pTj)3  ' 


(5.  13) 


and  we  shall  not  concern  ourselves,  in  what  follows,  with  the  complementary 
expressions  for  (7^  and  Mg. 

INTERIOR  AND  EDGE  ZONE  SOLUTION  CONTRIBUTIONS 

We  accept  as  known  the  fact  that  the  solution  to  be  obtained  will,  for 
sufficiently  small  values  of  p,  be  composed  of  two  contributions,  one  of  them 
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being  an  "interior"  contribution  g^,  F^,  with  characteristic  length  b = b.  = (ac)‘ 

and  the  other  an  "edge  zone"  contribution  g , F , with  characteristic  length 

e e ® 

b = b = c.  We  note  that  in  this  way  F.  = F.(4.)  and  F = F (£  ),  etc.  , but 
e 1 i i e e e 

we  shall  refrain  from  making  the  distinction  between  £.  and  £ explicit,  as  this 

i e 

will  cause  no  difficulty  as  the  analysis  proceeds.  With  £.,  £ and  T]  defined  in 

i e 

this  manner,  and  in  view  of  the  form  of  the  differential  equations  and  boundary 
conditions  for  and  F,  we  are  in  a position  to  stipulate  the  basic  order  of 
magnitude  relations 


g,.  g,.  F.,  F^,  g:.  g;,  g;.  g;.  F/.  f;  = 0(1)  , 


(6.1) 


and  we  shall  write,  in  place  of  equations  (5.2),  so  as  to  make  the  distinction 
between  the  two  solution  contributions  explicit 


(1  - = E('I'.g.  + '&  g ) , V = V.F.  + V F 

ii  60  li  cc^ 


(6.2) 


with  the  form  of  the  two  differential  equations  for  the  set  g.,  F.  and  for  the 

set  g , F being  distinct,  to  the  extent  that  this  is  required  to  account  for  the 
e e 

difference  between  b.  and  b . 

1 e 

Setting  b = b^  = (ac)^,  we  obtain  from  equations  (5.4)  and  (5.  5) 


Setting  b = b = c,  we  obtain  instead 
e 


(6.3) 


(6.4) 


J - 0 . 


p 1 - 0 e' 


'I'  [g"  - g’  + a^g"]  + V [(1  + PT7)F"  + F ] = 0 

e e 1 + prj^e  ^e-*  e*-  ' e 1 + pT)  e^ 


(6.5) 


(6.6) 


In  view  of  the  difference  between  £.  and  £ , the  functions  g.  and  g must 

I e i e 
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individually  satisfy  the  face  boundary  conditions 


rj  = ± 1:  g.  = 0 , gg  = 0 ' (6.7) 

and  for  the  same  reason  the  edge  conditions  at  infinity  (5.7)  must  be  satisfied 
individually, 

4 = «:  g;  = 0 , g = 0 , F.  = 0 , F = 0 . (6.8) 

I 6 1 e 


With  g.  and  F.,  and  g and  F , sofar  being  subject  to  a system  of 
11  e e 

uncoupled  requirements,  there  remains  only  the  system  edge  loading  conditions 
(5.8)  and  (5.9)  to  accomplish  the  necessary  coupling  for  the  determination  of  the 

i 

two  sets  of  functions.  We  find  that  equations  (5.8),  again  with  b^  = (ac)^  and 
b = c,  takes  on  the  form 


+'I'^g'(0,rj)]  + V.p^F.'(O)  + V^F^(O)  = -cu 


(6.9) 


and  equations  (5.9)  take  on  the  form 


V.F.(O)  + V F (0)  =0  , 

11  e e 

or 

(.1  1 

['I'P^g'(0,T7)  + •I'g'(0,T7)]dr)  = 0 . 
J 1 1 e e 

-1 


(6.  10a) 


(6.  10b) 


Having  the  relations  (6.2)  to  (6.  10)  our  next  step  is  an  appropriate 

disposition  of  the  four  scale  factors  'I'^,  'I'^,  V^,  V^.  We  begin  by  observing 

that  in  order  to  be  able  to  satisfy  the  face  boundary  conditions  (6.  7)  in  such 

a way  that  not  only  the  functions  g.  and  g , but  also  the  functions  F.  and  F 

1 e 1 e 

are  involved  it  will  be  necessary  to  have  the  terms  with  g * in  equations  (6.4) 
and  (6.6)  of  the  same  order  of  magnitude  as  the  terms  with  F"  in  these  equations. 
Accordingly,  we  get 

>1^.  = PV.  , 'I'  = V . (6.11) 

1 1 e e 

Introduction  of  (6.  11)  into  the  edge  conditions  (6.9)  and  (6.  10)  then  gives 


further 


with  (6.  10a)  remaining  unchanged,  and  equation  (6.  10b)  assuming  the  form 


pi  3 

J [V.P=gj^{0.T))  + V^g^(0,ri)]dTl  = 0 . 


(6.  13) 


Having  equation  (6.  12),  we  now  impose  the  requirement  that  both  interior 
and  edge-zone  solution  contributions  participate  in  the  satisfaction  of  this  non- 
homogeneous  edge  condition.  This  means  that  we  must  have 


V.P^  = V . 
1 e 


(6. 14) 


Finally,  we  introduce  a weighted  average,  0,  of  the  variation  of  edge 
rotations  in  thickness  direction  by  writing 


u = -0s(77)  , 
> * 


(6. 15) 


J (1  - r?®)s(T?)dT)  = j , 0 = -IJ  (1  - T)^)G  ^d7)  , 


(6.  16) 


(so  that  0=0  when  u = -0  = const. ). 

o , r o 

With  (6.  14)  and  (6.  15)  and  with  the  further  stipulation  that 


V,  = cp  , 


(6.  17) 


the  nonhomogeneous  edge  condition  appears  now  in  the  form 


(6.  18) 


and  the  homogeneous  edge  conditions  (6.  13)  and  (6.  10a)  become 


F (0)  + P^F  (0)  = 0 , 
1 e 


(6.  19a) 


J [pg/(0,7))  + g^{0,T7)]dr7  = 0 . 


(6.  19b) 


I 


Equations  (6.  18)  and  (6.  19)  are  to  be  used  in  conjunction  with  the 
consequences  of  introducting  (6.  11)  into  equations  (6.3)  to  (6.6),  that  is 
in  conjunction  with  two  systems  which  may  be  written  in  the  form 
1 


1-  1 ^2  1 X n 

J g.dn  + —^^- — -F 
-1 


1 - p 


, r ^ 

i ^ ^(1  + P7J)* 


= 0 


g + F + p 

“i  i 


and 


g:  - i^F. 

+ VF':  - -7 - 

1 1 + PTJ 


(6.20) 


= 0 


. , Vg  dn 

J ^ J imsi  ^ 


1 - 1 + p 


-in 


®e  " ®e  e ^ 


g. 


1 - p e| 
PUF 


1 + PV 


+ r)F  + 


e 1 + P7] 


= 0 


= 0 


(6.21) 


The  solution  of  (6.20)  and  (6.21),  subject  to  the  edge  conditions  (6.  18), 
(6.  19)  and  (6.8)  and  subject  to  the  face  boundary  conditions  (6.7)  is  to  be 
introdx'ced  into  equations  (6.2)  for  and  V.  Upon  taking  account  of  the  relations 
(6.  11),  (6.  14)  and  (6.  17)  we  have  then 


(1  - = EcP(g  + p2g.)  , V = v^8(F.  + P^F  ) . 


(6.22) 


Furthermore,  introduction  of  'P  and  V from  equation  (6.2)  into  equations  (5.  10) 
to  (5.  12)  gives  as  expressions  for  stresses 


e3 


g;  + p^g: 


a = 

X 1 - t'S  1 + PT) 


T = - 


E8 


K + 


1 - 1^2  1 + pry 


E0 


> J. 

g*  + p2g:  p2F.  + pF 

e i 1 e 


1^2  1 + prj 


1 + P77 


(6.23) 


(6.24) 


and  as  expressions  for  the  stress  couple  M and  the  stress  resultant  Q 

EBc’  r' , T 

“*  = -TTFj  (8^  + . 

-1 


E8c 
1 - u- 


J (?  ' + Pg.')d77  , 

e 1 


(6.25) 


(6.26) 


-1 
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We  note  that  while  the  above  formulas  for  M and  Q give  the 

X ® 

impression  that  the  edge-zone  solution  contribution  makes  a dominant 

contribution  to  these  quantities — which,  if  true  would  detract  much  from  the 

significance  of  conventional  two-dimensional  shell  theory — we  will  find  shortly 

that  J g d77  comes  out  to  be  of  a smaller  order  of  magnitude  in  p than  g.drj, 

“1  6 “1  1 

so  that  no  problem  of  this  sort  does,  in  fact,  arise. 

Beyond  this  we  may  also  note  that  the  differential  equation  for  g in 
(6.21)  is  effectively  of  second  order,  in  place  of  a fourth-order  bi-harmonic 
problem  which  is  encountered  in  the  analysis  of  the  corresponding  problem 
for  an  isotropic  homogeneous  shell  medium.  It  is  this  difference  in  the  order 
of  the  differential -equation  problem  for  the  edge-zone  solution  contribution 
which  makes  possible  an  exact  solution  of  the  complete  problem  for  the  case 
of  the  transversely  inextensional  medium. 

PARAMETRIC  EXPANSIONS  FOR  INTERIOR  AND 
EDGE-ZONE  SOLUTION  CONTRIBUTIONS 

As  inspection  of  the  different’^!  equations  (6.20)  and  (6.21),  in  conjunction 
with  the  boundary  conditions  (6.7),  (6.8),  (6.  18)  and  (6.  19)  indicates  the  possibility 

i f 

of  an  expansion  of  the  solutions  in  powers  of  the  parameter  p^. 


F = F + p^F  + pF  + . . . , g = + P*g  + Pg  + . . . , 

O 1 « 0 13 


(7.  1) 


and  observing  that  (llp)!!m[{l  + P)/(l  - P)]  = 2[l  + p^/3  + . . .],  we  have  then  from 
equations  (6.  20), 


g'  dTJ  + 2(1  - i/*)F  = 0 

J in  in 


» + F =0 

’in  in 


n = 0,1 


(7.2) 


/gr„d„t2(l  . -Fl'g  d.,  ; n = 2,  3 

'1  -1 


g!  + F.  = g:  - a^gf  - tjf!  - VF. 

in  in  in-a  in-a  in-a  in-a  ’ 


(7.3) 


Except  for  the  form  of  (6.  19a)  this  expansion  would  be  in  powers  of  P,  rather 

than  P^.  This  means  that  when  (6.  19b)  applies,  one-half  the  terms  in  our 
expansion  will  vanish  automatically. 
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with  analogous  relations  for  n = 4,  5,  . . . , which  will  not  be  utilized  in  what 
follows,  and  equations  (6.21)  give  in  the  same  manner 


.1 

I g'  dT)  = 0 . g"  + + F"  = 0 ; n = 0,  1 , 

J ®en  en  en  en 

-1 


(7.4) 


r^g*  dT?  = 0 . g"  +ar*g'  + F"  = g-  - T?F"  ; n = 2,  3 . (7.5) 

J *en  en  en  en  en-a  en-a 

-1 

[ g"  dT?  = r g dT?  - 2(1  - U^)F  . - 

J ®en  J en-4  en-*  ; n = 4,  5 . 


(7.6) 


The  associated  boundary  conditions  which  follow  from  (6.7)  and  (6.8)  are 
n . i 1:  8i„  = 0 • *.n  = “ ■ 

« = *i„  = 0 • = 0 • (’-81 


The  coupling  between  interior  and  edge  zone  contributions  will  be  effected 
in  the  process  of  satisfying  the  boundary  conditions  for  the  loaded  edge  of  the 
shell  which  follow  from  (6.  18)  and  (6.  19)  in  the  form 

F/  (0)  + F'  (0)  + (0,T?)  = s(77)  . 


lO 


eo 


eo 


f/  (0)  + F'  (0)  + a^g'  (0,T?)  = 0 , 

11  Cl  Cl 


(7.9a) 

(7.9b) 


^iV°^  ^en^°^  ^ " “ ^®en-a^° ' “ ^in-a^® ' " 2,  3.  (7.9c) 


and 

F.  (0)  = 0 , F.  (0)  + F (0)  = 0 ; n = 1,  2,  3,  . . . 
lo  in  en-i 


or 

..1 


(7.  10) 


I g'  (0,T?)dT?  = 0 ; n = 0.  1 , 

J en 
-1 


(7.  11a) 


1 


It  is  in  this  one  place  that  we  require  equations  for  n = 4,  5,  in  order  to  complete 
the  determination  of  our  expansions  terms  upto  n = 3. 
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l-l 

J ■ 'J  ' " = 2,  3,... 

-i  -1 


(7.  11b) 


DETERMINATION  OF  INTERIOR  SOLUTION  CONTRIBUTION  TERMS 

We  begin  by  deternnining  the  interior  solution  contribution  to  the  extent 
that  this  is  possible  without  bringing  in  the  fonn  of  the  associated  edge  zone 
solution  contribution.  In  doing  this  we  will  arrive  at  a set  of  ordinary  differen- 
tial equations  for  the  dimensionless  deflection  functions  F.  . with  F beine 

in  lo  “ 

equivalent  to  the  corresponding  function  in  accordance  with  conventional  two- 
dimensional  shell  theory,  with  F^  accounting'for  the  effect  of  deviating  from 
the  "conventional”  assumption  of  an  T7-independent  edge  slope  function  s(77),  and 
with  F^  accounting  for  the  effects  of  transverse  shear  deformation  and  of 
element  width  changes  in  the  same  manner  as  previously  found  in  [l]. 

We  complement  these  results  by  the  discovery  of  the  complete  system 
of  boundary  conditions  for  for  an  arbitrarily  prescribed  edge  slope  function 

s(T)),  through  use  of  the  differential  equations  for  the  edge-zone  solution  contri- 
butions, ^’ut  without  having  to  determine  any  part  of  the  edge -zone  portion  of  the 
solution  of  the  complete  problem. 

Considering  the  fact  that  F^^  = F.^(C),  we  find  from  equations  (7.2),  in 


in 


conjunction  with  the  appropriate  statement  in  (7.7)  the  simple  result 


g; 


- ^ 


in 


in 


IV  a 

F:  + 3(1  - = 0 

in  m 


n = 0,  1 


(8.  1) 


In  the  same  way  we  find  from  equations  (7.  3)  and  (7.  7),  and  with 


(1  - l/V®  = E/G, 


- V _ E 5 - 677 


g,„  = rr  1 3^  f: 

in  2 in  3 in-a  \ 

C * - ‘'’If'in  = (It  - > " = 2.  3 . 


1^)f. 

/ in- 


(8.2) 


We  now  turn  to  the  form  of  the  boundary  conditions  for  the  functions  F. 

We  find  from  (7.8)  the  obvious  requirements  of  vanishing  F (and  F^  ) for 

in  in 

5 = ».  As  far  as  conditions  for  ^ =0  are  concerned,  we  have  as  one  of  the  two 
conditions  for  F.^,  for  one  of  the  two  cases  considered  in  (7.  10)  and  (7.  11),  that 
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F.  (0)  =0.  The  other  conditions,  in  accordance  with  (7.9)  to  (7.11)  all  appear 

lO 

to  require  the  simultaneous  determination  of  F and  g . We  will  show  next  how 


en 


en 


to  avoid  this  complication  to  a significant  extent. 


DERIVATION  OF  BOUNDARY  CONDITIONS  FOR  F.  , WITH  OR  WITHOUT 

in 

INVOLVEMENT  OF  EDGE -ZONE  SOLUTION  CONTRIBUTION 

We  begin  by  observing  that  equations  (7.4)  and  (7. 5),  in  conjunction  with 
the  conditions  g (“,^)  = =0  imply  the  relations 


en 


en 


J = J g;„(?.Wd,7  = 0 , 


for  n = 0,  1,  2,  3.  Furthermore,  from  (6.21), 

1 ® 

I g^(«.ri)dil  = (I  - f F^^(?)de  . 


(9.  la) 


(9.  lb) 


-1 


'.V’th  the  help  of  (9.  1),  and  through  use  of  the  first  relation  in  (8.  1)  and 
(8.2),  we  may  write  the  set  of  conditions  in  (7.  10)  and  (7.  11)  in  the  form 


o 

o 

= 0 , 

or 

F-fO) 

= 0 , 

F^(0) 

= -F  (0)  ' 

eo 

or 

= 0 , 

= 

en-1 

or 

/6E 
= \5G 

(9.2a) 

(9.2b) 


in 


(0) 


- 3(1  - y®)  r F (C)d4 

J en-2 


n = 2,  3 


(9.2c) 


It  remains  to  establish  a second  condition  at  ^ = 0 for  F.  , through  use 


in 


of  equations  (7.9),  with  a view  towards  letting  this  condition  to  be  as  much  as 
possible  a condition  for  terms  of  given  quantities.  In  order  to  accom- 

plish this  purpose,  we  consider  an  integrated  version  of  the  edge-zone  differential 
equations  (7.4)  and  (7.5),  with  an  appropriate  weighting  function,  as  follows. 

We  write,  as  a consequence  of  equation  (7.4), 
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n 


0,  1 . 


(9.4) 


J'  (I  - t = 0 , 

-1 


In  this  we  write  first. 


r (1  - T7’)g”  dr?  = [{1  - T]=)g-  + 2T)g  ]'  - 2 1 g dr)  . 

J en  en  j en 

-1  -1  -1 


j 

■1  -1 

with  the  conclusion,  which  follows  from  (9.  1)  and  (7.8),  that 


f = 0 • 

•I  en 

-1 


(9.5) 


Introduction  of  (9.  5)  into  (9.4)  and  observation  of  the  edge  conditions  at  infinity- 
then  gives  further 


r 


-1 


(1  - T)  )(«^g'  + F'  )dT7  = 0 ; n = 0,  1 . 

en  en 


(9.6) 


Upon  using  (9.6),  in  conjunction  with  (7.9a,b),  we  obtain  as  conditions  for 


F.  and  F.  , 
lO  u 


f'  (0)  = 1 , F'  fO) 

LO  U 


(9.7a,  b) 


and  we  see  that,  in  fact,  the  determination  of  F.  may  in  all  cases  considered 

LO  

here  be  carried  out  without  explicit  reference  to  the  edge  zone  contribution  to 
the  complete  solution  of  the  problem. 

It  suggests  itself  to  see  to  what  extent  application  of  the  same  procedure 

to  the  differential  equations  (7.5)  may  make  possible  a reduction  of  the  boundary 

conditions  (7.9c)  which  involve  F^^  and  F.\  We  find  first,  proceeding  in  the 

same  way  as  in  going  from  (9.4)  to  (9.6),  and  observing  that  rjF"  is  an  odd 

en-a 

function  of  T) , that  now 


J 


1 

(1  - r)®)(Qf^g''  + F"  )dT)  = 2 r)g  dT)  , 

en  en  J en-a 

-1  -1 


(9.8) 


for  n = 2,  3. 
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Evidently,  equation  (9.8),  together  with  the  edge  conditions  at  infinity, 
allows  us  to  conclude  that 


J (1  - T?^)[«®ggjj(0.T7)  + F^^{0)]dT)  = -2  j T][ 


-1 


-1  o 


(t77)dCdr7  , 

en-a 


(9.9) 


for  n = 2,  3.  Therewith,  and  with  the  consequences  of  the  first  relation 
in  (8.  1),  we  obtain  altogether  that 


3 ^ QB 

-1  o 

3a®  3 / 

+ — J (n-T)=’)g^^_^(O.T,)d77  , 


(9.  10) 


for  n = 2,  3.  We  note  as  a very  important  consequence  of  (9.  10)  the  fact  that 

when  g is  an  even  function  of  77,  which  is  the  case  when  s(T7)  is  even,  then  the 
^eo 

integral  terms  on  the  right  vanish  for  n = 2 and  we  have  a boundary  condition  for 
F.  as  well  which  is  free  of  any  reference  to  the  edge-zone  solution  contribution. 


— la 


EXPANSIONS  FOR  DEFLECTION,  AXIAL  BENDING  MOMENT  ^ND 
TRANSVERSE  SHEAR  STRESS  RESULTANT 

In  order  to  see  the  possibility  of  determining  M and  Q,  and  also  V,  for 


the  conventional  displacement  boundary  condition  case  of  an  TJ- independent  u 
upto  terms  of  relative  order  p,  solely  on  the  basis  of  a determination  of  the 
interior  solution  contribution  we  rewrite  our  earlier  expressions  for  these 
quantities  as  follows. 

On  the  basis  of  equations  (6.22''  and  (7.  1), 


, r 


11  1 

V = -8c2a^[F.  + P^(F.  +F  ) + p(F.  +F  ) + ...]  . 

lo  11  eo  U2  ei 


(10. 1) 


On  the  basis  of  (6.25),  (6.26),  (7.  1),  in  conjunction  with  equations  (8.  I] 
and  (8,2)  for  the  g^^  and  equations  (9.  1)  for  the  g^^^. 


M 


2 

■3 


E3( 


(1  - I^“)a" 


r{Fr  + p^F.''  + P 
* 10  11 


(10.2) 
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2 E3c^  1 

■3  (1  - l/“)a  I 


F"'  + P^F"' 
10  u 


We  note  that  it  will  be  possible  to  determine  V,  and  Q.  upto  terms  of 

relative  order  P,  soley  on  the  basis  of  a consideration  of  the  interior  solution 

contributions  F . F and  F . to  the  extent  that  it  is  possible  to  determine  these 
10  ii  is 

functions  without  explicit  consideration  of  the  edge-zone  solution  portion,  as 

discus  d on  page  21,  and  to  the  extent  that  F and  F turn  out  to  vanish 

GO  d 

altogethe  r . 

The  same  conclusions  apply  insofar  as  the  determination  of  the  stiffness 
coefficients  and  is  concerned,  inasmuch  as  we  have  that 


Kpp  = M^(0)/3  , Kyg  = -Q(0)/3 


(10.4) 


when  V(0)  = 0,  and 

Kg  = M^(0)/S  , 


when  Q(0)  = 0. 


(10.  5) 


EXPANSIONS  FOR  STRESSES 

Introduction  of  the  series  expansions  (7.  1)  into  equations  (6.23)  and 
(6.24)  gives  as  expansions  for  the  stresses  T,  Og,  upto  terms  of  relative 
order  P, 


+ p"(g:,  + g-L)  + p(g;3  + g'ix  - ^ - ^ ' 


X \ ^®eo 


= - — ~~ — fg ' + + Pis'  + g-^  ~ 'Hs'  )+•••]  > 

1 - patSeo  ®ei  '^ea  ®io  '®eo' 


ei  10 


On  = VO  + E2[p^F.  + P(F.  +F  ) + ...]  , 

0 X 10  u eo' 


(11.  1) 


(11.2) 


(11.3) 


^ Note  that  (10.4)  and  (10.5)  are  consistent  with  (4.1)  and  (4.2)  in  the  light  of  the 
defining  relation  (6.  I6)  for  3. 
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Jl 


r^. 


4 


1 


We  note  that,  in  contrast  to  the  results  for  deflection,  axial  moment 

and  transverse  resultant,  the  distribution  of  stress  in  the  shell  will,  for  the 

general  case,  be  such  that  edge  zone  contributions  come  out  to  be  of  a higher 

order  of  magnitude  than  interior  region  contributions.  This  is  not  so  for  the 

"exceptional"  cases  for  which  g^^  and  g^^  vanish  identically.  The  significance 

of  this  distinction  will  become  apparent  in  our  discussion  of  the  exceptional 

case  of  an  TJ- independent  u , in  comparison  with  the  case  where  u is  assumed 

* ^ , r 

to  be  a quadratic  function  of  TJ. 

SOLUTION  OF  INTERIOR  AND  EDGE- ZONE 
DIFFERENTIAL  EQUATIONS 


Determination  of  the  interior  solution  contribution,  in  accordance  with 
equations  (8.  1)  and  (8.2),  is  the  same  as  in  the  earlier  work  concerned  with  this 
contribution  only  [l].  We  find,  upon  taking  account  once  for  all  of  conditions 
for  4 with  4m*  = 3(1  - U^),  that 


F. 

in 

= (A  cosm^  + B sinm4)e  ^ ; n = 0,  1 

n n 

f 

(12.  1) 

and 

F. 

in 

-m^ 

= (A  cosmf  + B sinm4)e  ’ 
n ^ n 

■ (sG  ■ ' n - 2,  3 , 

(12.2) 

with  the  associated  functions  g^^  following  from  the  above  in  accordance  with 
(8.  1)  and  (8.  2). 


Having  F.  and  g.  , as  indicated  above  we  now  turn  to  the  determination 
in  in 

of  F and  g , for  n = 0,  1,  2,  3,  in  accordance  with  equations  (7.4)  to  (7.8). 
en  ^en  \ i 

Writing  as  before  as  differential  equations  for  g and  F , the  set 

en  en 


g"  + + F"  = 0 ; n = 0,  1 , 

®en  ®en  en  ' 


g"  +«"g''  +F"  = g*  -TJF"  ; n = 2,  3 , 

en  en  en  en-a  en-a 


(12.3) 

(12.4) 


we  take  account  of  the  first  relations  in  (7.4)  and  (7.  5),  in  conjunction  with  the  g 

conditions  for  ^ = ®,  so  as  to  establish  as  a system  of  three  boundary  conditions  9 
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for  rj  = ± 1 , 


J g„(S,Wd^  = 0 , g,„(C,±i)  = 0 . 


(12.5) 


with  the  presence  of  the  terms  (4)  in  (12,3)  and  (12.4)  making  it  possible  to 

en 

satisfy  the  three  conditions  in  (12.5),  when  without  these  terms  the  first  condi- 
tion (12.5)  would  rule  out  significant  portions  of  the  solutions  of  (12.3)  and  (12.4). 


To  see  that  this  is  in  fact  so  we  write  (12.3)  and  (12.4)  in  the  differentiated 


form 


(12.6) 


•••  3 //  • ••  A 

g+«g  =g  -F  ;n=2,  3, 
®en  en  en-3  en-s 


(12.7) 


and  determine  first  the  solution  of  (12.  6). 


We  find,  by  separation  of  variables,  that  equation  (12.6)  has  particular 


solutions 


g = (C  + C cos  Xrj)e  +(Dsmur))e 

*enp  'la" 


(12.8) 


where  C^,  C^,  D,  X and  U are  arbitrary  constants.  Satisfaction  of  the  three 
boundary  conditions  (12.5)  requires  that  these  arbitrary  constants  satisfy  the 


relations 


C + C X ^ sinX  = 0 , C,  + C cos  X = 0 , Dsinli  = 0 
la  la 


(12.9) 


The  conditions  for  non-trivial  solutions  in  (12.9)  are 


tan  X = X , sinli  = 0 , (12,10) 

for  all  positive  roots  X = X.  and  ^ = W,  for  j,k  = 1,  2,  3 It  is  evident  that 

1 ^ 

the  functions  sinu^77  form  an  orthonormal  system,  such  that  every  integrable 
odd  function  of  T]  in  the  interval  - 1 < T)  < 1 may  be  expanded  in  terms  of  them. 
Additionally,  it  can  be  shown  that  the  set 


cos  X.  - cos  X.T) 


r 


also  forms  an  orthonormal  system,  in  such  a way  that  every  integrable  even 
function  f(t7),  in  -1  < TJ  < 1,  with  the  additional  property  J fd77  = 0,  can  be  expanded 


in  terms  of  them. 


We  use  the  above  properties  to  construct  the  series  solution 


CD 

g = ^ C.  (fi.(V)e 


for  n = 0,  1. 


-X.C/ff 


+ D,  (s  in  a r])e 
kn  k 


(12.  12) 


Before  proceeding  to  the  analogous  solution  of  (12.7)  we  return  to  equation 
(12.3)  in  order  to  establish,  with  the  help  of  the  conditions  F^^(oo)  = («)  = 0,  as 

solution  F (^)  which  is  associated  with  g in  (12.  12), 


en 


en 


en 


en 


en 


; n = 0,  1 . 


(12. 13) 


We  now  consider  equation  (12.4)  where  it  remains  to  determine  the 
appropriate  particular  solution  corresponding  to  a right-hand  side  which  can 
be  seen,  with  the  help  of  (12.  11)  to  (12.  13),  to  be  of  the  form 


/sinX.T?  \ 

y^C.  - f-  + 

A-/  jn-3 \ cos  A.  j 7 


-X.4/a 


+ D,  (U,  cos  u T))e 
kn-3  k k 


We  find  that  a particular  solution  of  (12.4)  with  the  above  expression  in  place  of 
the  original  right-hand  side,  which  also  satisfies  the  three  boundary  conditions 
(12,5),  is  such  that  altogether 


-xe/“ 

g = y^C.  cp.(T7)e  ■’  + D,  (sinU,T7)e 

®en  ^ jn  j kn'  k ' 


T'c 

jn- 


rjcosX.T)  sinX.T)  \ -u  f/of 

jL  . i — . ] ) g ^ 

3\X.  2X.cosX.  2X.sinX./ 

J J J J j/ 


We  note  that  it  is  this  orthogonality  property  of  the  functions  (p,  which  is  responsible 

for  the  advantages  of  the  present  solution  of  the  general  problem,  in  comparison 
with  the  solution  of  the  cor respondiag  problem  for  the  isotropic  shell. 
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TJ  sinU.  T?  - 

K 


COS  Uj^77  - cosUj^\  -u^4/“ 


(12.  14) 


-X.^/or 


EX.  ''C  . e ^ ^ ' u ^ cos  U.  D e 

j jn  2 k K kn-3 


-U  C/Qf 

k 


for  n = 2,  3. 


(12. 15) 


It  now  remains  to  determine  the  constants  of  integration  C . and  D, 

jn  kn 

in  conjunction  with,  or  preferrably  subsequent  to,  the  determination  of  the  constants 

A and  B in  the  interior  solution  portion  F.  as  given  by  (12.  1)  and  (12.2). 
n n in 

DERIVATION  OF  BOUNDARY  CONDITIONS 

FOR  DETERMINATION  OF  CONSTANTS  IN  g 

en 

We  begin  by  recalling  that  the  complete  system  (7.9)  to  (7.  11)  of 

conditions  at  the  edge  C = f . = 0 of  the  shell  had  be*;.!  reduced,  insofar  as 

possible,  to  a system  of  conditions  for  the  F,  alone,  consisting  of  equations 

(9.2),  (9.3),  (9.7)  and  (9.10).  A characteristic  of  these  conditions  is  that  in 

them  the  effect  of  the  edge-zone  solution  contribution  manifests  itself  through 

the  presence  of  certain  definite  integrals  with  respect  to  ry  of  g (4»^)»  but  with- 

en 

out  appearance  of  the  g^^  themselves.  Consequently,  these  conditions  do  not 

contain  those  elements  of  the  complete  system  of  edge  conditions  (7.9)  to  (7.  11) 

which  are  required  for  a determination  of  the  series  coefficients  C and  D in 

jn  kn 

Sen- 

A re-inspection  of  equations  (7.9)  to  (7.  11)  makes  it  clear  that  the 
remaining  supplementation  of  the  boundary  conditions  (9.2),  (9.3),  (9.7)  and 
(9.10)  must  come  from  equations  (7.9),  and  not  from  (7.  10)  and  (7.11). 

We  obtain  the  first  of  the  remaining  transformed  conditions  by  introducing 

(9.7a)  into  (7.9a)  and  by  eliminating  F^  (0)  from  this  relation  by  considering  that 

eo 

J*  g'  dr)  = 0.  This  gives 
-1  eo 


(13.  1) 
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We  finally  introduce  f/  (0)  from  (9.  10)  into  (7.  9c)  and  consider  that 
1 1 

here  J*  g'  drj  = o,  and  g/  (0,TJ)  =-7(1  - T7®)F"  (0).  Therewith  equation  (7.  9c) 

-1  en  in-3  Z m-a 

becomes 


- I T’ 

-1 


T7g_  (O.n)dv  - (0) 

en-3  Z 3 in-3 


(13.3) 


for  n = 2,  3. 

Having  equations  (13.  1)  to  (13.  3),  in  conjunction  with  equations  (12.  12) 
and  (12.  13)  for  g and  F , we  are  now  in  a position  to  arrive  at  some  general 

6X1  6n 

conclusions,  without  additional  analysis.  To  whit 

(1)  The  functions  g and  F will  vanish  identically,  for  all  edge 

ei  ei 

conditions  cases  here  under  consideration. 

(2)  The  functions  g and  F will  vanish  identically  for  the  special  case 

eo  eo 

s(T7)  = 1,  but  will  otherwide  have  a dominant  effect  in  equations  (11.1)  to  (11.  3) 

j. 

for  the  distribution  of  stress,  and  a direct  effect  of  relative  order  on  the 

values  of  the  deflection  V,  in  accordance  with  equation  (10.  1),  as  well  as  an 

1. 

indirect  effect  of  relative  order  p^  on  and  Q,  in  accordance  with  equations 
(10.2),  (10.  3)  and  (9.2b). 

(3)  The  functions  g and  F will  make  a contribution  of  relative  order 

es  es 

p on  stresses,  in  accordance  with  (11.  1)  to  (11.  3)  but  will  make  no  such  contri- 
bution to  V,  M and  Q. 

X 

THE  CASE  OF  A UNIFORM  EDGE  ROTATION 


We  assume  now,  as  in  [4],  an  axial-edge  displacement  u = - a-), 

giving  u = -3  and,  in  accordance  with  (6.  16),  3 = -(3/4)(-3  )(4/3)  = 3 and 
, r o 00 

therewith,  in  accordance  with  (6.  15), 
s = 1 , 


(14.  1) 


and  then,  as  discussed  in  the  preceding  section 


g = 0 , 
eo 


F =0 
eo 


(14.2) 


With  (14.  2),  and  with  g =0andF  =0,  we  have,  from  (9.2),  (9.3),  (9.7) 

ei  ei 

and  (9.  10)  as  edge  conditions  for  the  interior  solution  contribution 
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F*'(0)  = F"'(0)  = F:‘'(0)  = F"'(0)  = 0 

lO  U 13  13 


f'  (0)  = 1 . F'  fO)  = F'  iO)  = f'  (0)  . (14.3) 

lO  11  13  13 

together  with 

= ° ' 

for  the  case  V(0)  = 0,  and 

Fr^(O)  = F|^(0)  = F'^(O)  = F"'(0)  = 0 (14.5) 

for  the  case  Q(0)  = 0,  with  equation  (14.  3),  together  with  (14.  4)  or  (14.  5),  to  be 

used  for  the  determination  of  the  constants  of  integration  A and  B in  (12.  1)  and 

n n 

(12.2). 

We  omit  the  elementary  calculations  leading  to  the  appropriate  values 

of  the  A and  B . With  their  help  there  follows 
n n 

i 1 

V(0)  = -c^3i^?[F.  (0)  + PF.  (0)  + . . .] 

lO  13 


M (0)  = - 

X 


2m  a 
Ec^ 


m“  \20O  4/  J ' 

5 

^3|f;'  (0)  + p[f;'(0)  - (0)1  + ..  J 

* 2 ( lO  1 \5G  / lo  J I 


(14.6) 


Ec^  af,  0 /3E  u\ 


(14.7) 


for  the  case  Q(0)  = 0,  and 


M (0)  = -^3  1 - 

X'  ' 3 ^ 

m a^ 


■°  1 

/9E 

m^' 

^20G 

+ "T  I + • • • 


for  the  case  V(0)  = 0 where,  it  is  recalled,  4m*  = 3(1  - and  P = c/a. 


(14.8) 


(14.9) 


We  note,  specifically,  that  the  results  in  (14.6)  to  (14.9)  have  been  obtained 
without  determination  of  any  edge*zone  solution  contribution  and  that  furthermore 


the  stiffness  coefficients  and  which  are  associated  with  (14.  8),  in 

accordance  with  (10.4),  as  well  as  the  coefficient  Kq  associated  with  (14.7) 

P 

in  accordance  with  (10.5)  are  consistent  w ith  our  earlier  bound  results  in  [4], 
upon  specializing  these  so  as  to  correspond  to  the  limiting- type  assumption  of  a 
medium  unable  to  experience  transverse  normal  strains. 

Having  determined  the  we  may  now  obtain  the  here  leading  terms 

g^  and  of  the  edge -zone  solution  portion  of  the  complete  solution  of  the 
problem  from  equations  (12.  14)  and  (12.  15),  with  = 0 and  D.^  = 0,  in 
conjunction  with  the  boundary  condition  (13.3)  which  now  reduces  to  the  form 

(0)  , (14.10) 

r lo 

or,  equivalently,  to 


gJO.V) 

€3 


Hi  ■ 


io' 


Vx.c.  ©.  + u.  D,  Sinu,  T)  = - Ijf:"  (0) 

J J3  J k ka  k 2 \ 3/  lo  ' 


(14.  11) 


in  the  interval  - 1 < Tj  < 1.  The  orthonormality  properties  of  the  functions  and 
sinUj^T)  then  give 


= 0 , Xj 


ka 

and  therewith. 


,.C.  = jF:'(0)f  (t)^  - l-j'p.dT]  = |^f:'(0) 

J ja  2 lo  ' 3/  J 3Xj  lo 


(14. 12) 


8 = 

'’ea  3 J 


(T?)e 


-X.^/Of 
J 


(14.  13) 


A consideration  of  (14.5)  now  gives  that  g vanishes  throughout,  just  as 

ea 

g and  g , for  the  case  Q(0)  = 0.  For  the  case  V(0)  = 0 the  function  F.  , with 
eo  ei  lo 

F.^  (0)  = 1 and  F.  (0)  = 0,  gives  the  relation  F^  (0)  = 2m'. 

lO  lO  lo 

Having  determined  g it  becomes  possible  to  evaluate  the  distribution  of 

ea 

stress  in  the  shell,  in  accordance  with  equations  (11.1)  to  (11.3),  with  these 
equations  reducing,  for  the  present  problem  to  the  form 


a 

X 


1 - 1/3 


+ 


T 


_£i£_rg' 

1 - i/3^^io 


+ 


(14.4) 
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with  a corresponding  expression  for  1 q.  We  see  once  again  that  the  "elementary" 

interior  contribution  in  (as  well  as  in  o q)  dominates  the  supplementary  edge 

zone  contribution.  At  the  same  time  the  edge  zone  contribution  is  of  the  same  order 

of  magnitude  as  the  interior  contribution  insofar  as  the  transverse  shearing  stress 

T is  concerned.  We  note,  in  particular,  the  possibility  of  writing,  on  the  basis 

(14.  14)  and  (14.  10),  as  expression  for  the  edge  shear  distribution  for  the  case  V(0)  = C 


^ E0O  m^ 
T{0.r7)  « y- 


(14.  15) 


with  the  second  term  inside  the  braces  representing  the  edge  zone  effect,  which 
has  the  expected  property  of  making  no  contribution  to  the  shear  stress  resultant 

Q. 

A CASE  OF  NON-UNIFORM  EDGE  ROTATION 

We  now  consider,  as  an  example  for  which  the  results  of  the  standard 
two-dimensional  shell  theory  are  complemented  in  an  essential  way  by  three- 
dimensional  considerations,  the  case  for  which  the  axial  edge  displacement  is 
prescribed  in  the  form 


u = -0 


f 1 k (r  - a)^ 
(r  - a)  + J 


u ^ = -Ml  + kT)=] 
, r o 


(15. 1) 


and  therewith,  in  accordance  with  (6.  16)  and  (6.  15), 


P = -P 


1 + kT)^ 
1 + k/5 


(15.2) 


^(l+|k)  , s(T,) 

We  can  now,  as  before,  determine  the  first  term  of  the  interior  solution 


contribution,  with  the  help  of  the  edge  conditions 

F'  (0)  = 1 , F.  (0)  = 0 or  F*' (0)  = 0 (15.3) 

io  to  — lo 

which  follow  from  (9.  2a)  and  (9.  7a) , depending  on  whethe r V(0)  =0o£Q(0)  =0 
is  prescribed,  and  we  note  that  this  determination  does  incorporate  information 
on  the  nature  of  the  shape  functions  in  (15.  1),  by  way  of  the  defining  relation  (15.2) 
for  P. 

Introduction  of  g from  equation  (12.  12)  into  the  relation 
eo 


-30- 


(15.4) 


«*g'  (0)  = ~ 

*eo'  ' 1 + k/5 


which  follows  from  (13.  1)  in  conjunction  with  (15.2)  now  gives 


C. 

JO 


k/3 


ofX^  1 + k/5 
J 


D.  = 0 
ko 


(15.5) 


Having  g , and  g.  , and  recalling  that  g^  = 0 for  all  cases,  we  now  have 


eo 


'lO 


as  expressions  for  stresses,  upto  and  including  terms  of  relative  order  p^, 


-E0  / 


1 - 1/2  ^eo  ' 


(15.6) 


as  well  as  CTg  = l/a^  + where,  notably,  all  components  of  stress 

are,  within  the  narrow  edge  zone,  of  a higher  order  of  magnitude  than  the 

exoected  stresses  associated  with  the  interior  solution  contributions  e.  and  F.  . 

lO  LO 

We  omit  listing  explicit  formulas  for  stresses,  and  instead  consider 

equations  (10.  1)  to  (10.  3)  for  V,  M and  Q,  where  we  note  that  evaluation  of  the 

i * 

contributions  of  relative  order  p^  and  P involves  the  relations 


TT  mL  -^k/3 

^eo^°^  ■ ■ 1 + k/5 
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I'  _ jf.  4Qt®k/3^-H  1 

J ^eo*^^  ~ * 1 + k/5ZF 

o J 


(15.7) 


Upon  evaluation  of  (10.  1)  to  (10.  3),  we  now  obtain  in  generalization  of 
equations  (14.6)  to  (14.9) 


V(0)  J_  j,  . Pj  V — 

^ca  ~ 2m)  m hc  l+k/5^X^ 


M (0) 
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1 

b1 

, /E_ 

2k  ^ 1 p 

/9E 

m 

''3G 

1 + k/5Z.^X"  ■ m=  ' 
J 

(20G 

1 

sl 

/H 

4k  y.  1 

m ' 

nc 

1 + k/5Z^X^ 

J 

P 

1 1 V M 1 

V 

Til 

5G  \ 

^ 1 + k/5  2^  Xj/ 

2 ^ 

(5.  JO) 


(15.  11) 


when  V(0)  = 0. 

In  order  to  evaluate  (15.8)  to  (15.  11),  we  note  that  with  the  successive 
roots  Xj  « 4.49,  7.72,  10.90,  14.07,  17.22,  20.37,  23.52,  26.67,  29.  81, 
32.99,  (j  + 7)ff,...  of  equation  (12.  10)  the  two  sums  of  negative  powers  of  X^ 
come  out  to  be 


^X"^  « 0.015  , « 


« 0.0029 


(15.  12) 


Insofar  as  the  interpretation  ^,5  equations  (15.8)  to  (15.  11)  is  concerned,  it 
should  be  noted  that  the  leading  terms  on  the  right  may  be  considered  as  equivalent 
to  the  consequences  of  ordinary  thin  shell  theory,  in  conjunction  with  the  solution 
of  the  problem  of  how  to  introduce  an  appropriate  representation  of  the  displace- 


ment condition  (15.  1)  into  this  theory.  Additionally,  we  find  that  while  the  effect 
of  transverse  shear  deformability  and  of  cross-sectional  width  changes  comes 
out,  as  expected,  to  be  of  relative  order  P,  the  effect  of  T7  -term  in  u ^ comes 
out  to  be  of  order  P*,  with  the  numerically  largest  values  of  these  correction 
terms  resulting  upon  letting  k tend  to  infinity,  with  a finite  limiting  value  of 

To  obtain  an  impression  of  the  numerical  consequences  of  replacing  a 

linear  distribution  u = u (r  - a)/c  by  a pure  cubic  distribution  u = u (r  - a)^/c^ 

o o 

we  set  in  equations  (15.8)  to  (15.11)  k = ® and  3 = -S^k/5  = 3u^/5c.  The rewith , 
and  with  equations  (15.  12)  there  follows 


2m  I m OG  m^  [20G 


(1  + 0.2)  +- 


(15. 13) 
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Ec^  Zm^U\ 


p 

3E 

m2 

20C 

(1 


0.2) 


when  Q(0)  = 0,  and 


(15.  14) 


(15.  15) 


(15.  16) 


when  V(0)  = 0. 

1 

As  might  be  expected,  the  shape  correction  terms  with  come  out  to  be 
numerically  quite  significant  for  moderately  thin  shells,  say  for  0 =0.1.  Addi- 
tionally, due  to  the  change  from  a linear  to  a cubic  edge  displacement  distribution 
a significant  modification  of  the  terms  with  P is  seen  to  occur  in  some  of  the  above 
expressions. 
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